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Essential Electromagnetism: Solutions 


Preface 


Preface 


This book gives the solutions to the exercises at the end of each chapter of my book “Essential 
Electromagnetism” (also published by Ventus). I recommend that you attempt a particular 
exercise after reading the relevant chapter, and before looking at the solutions published here. 
Often there is more than one way to solve a problem, and obviously one should use any valid 
method that gets the result with the least effort. Usually this means looking for symmetry in the 
problem - for example from the information given can we say that from symmetry arguments 
the field we need to derive can only be pointing in a certain direction. If so, we only need to 
calculate the component of the field in that direction, or we may be able to use Gauss’ law or 
Ampere’s law to enable us to write down the result. In some of these exercise solutions the 
simplest route to the solution is deliberately not taken in order to illustrate other methods of 
solving a problem, but in these cases the simpler method is pointed out. 

The solutions to the exercise problems for Each chapter of “Essential Electromagnetism” are 
presented here in the corresponding chapters of “Essential Electromagnetism - Solutions”. 

I hope you find these exercises useful. If you find typos or errors I would appreciate you 
letting me know. Suggestions for improvement are also welcome - please email them to me at 
protheroe.essentialphysics@gmail.com. 

Raymond J. Protheroe, January 2013 

School of Chemistry & Physics, The University of Adelaide, Australia 
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Electrostatics 


1 Electrostatics 


1-1 The surface of a non-conducting sphere of radius a centred on the origin has surface charge 
density cr(a, #,</>) = crocos# and is uniformly filled with charge of density po- Find the 
electric field at the origin. 

Solution 



At the centre of the sphere the electric field due to the volume charge will be zero because 
the contribution of a volume element located at r' will be exactly cancelled by that of an 
equivalent volume element at — r', so we only need to consider the surface charge. 


E(0,M) = 

47T£o J a z 


a* 

i r2ir n'K 

-~ / / cr(a, 0, 0)(—r) [a 2 sin 

'£ocl z Jo Jo 


^ /*27T r 1 

= - / dcf) d cos 9 (<jo cos 9) (—r). 

47T£o Jo J -1 


(1.1) 

( 1 . 2 ) 

(1.3) 


Because of the symmetry of the problem, the electric field at the centre can only be in the 
±z direction, and so we only need to find the z-component 


^ /*27T r 1 

E(0, 0,0) • z = -- / dcj) dcos9(<Jocos9)(-r)-z, 

47T£o Jo J -1 

1 f 1 

= - -27r / dcos9(ao cos9)(— cos0), 

4tt£:o 7_i 


(1.4) 

(1.5) 
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°o_ 

3^o 


(1.6) 


E(0 ,9,<f>) = --p-z. 


(1.7) 


1-2 A spherically symmetric charge distribution has the following charge density profile 


P(r, 0 , 4>) 


Po (' r < a) 

po(r/a)~P (r > a) 


(1.8) 


where /? is a constant (2 < f3 < 3). Find the electric field and electrostatic potential 
everywhere. 

Solution 

The charge density is spherically symmetric, with no dependence on 6 or 0, so the electric 
field must be in the radial direction and depend only on r. This is the ideal case to exploit 
Gauss’ law in integral form, 



(1.9) 


For r < a 



( 1 . 10 ) 


For r > a 



1 4 
£ 0 3 



( 1 . 11 ) 


1 4 3 

— oTTO p 0 + 


£ 0 8 


AirpoaP r(r') 3 ^"| r 


( 1 . 12 ) 



(1.13) 



(1.14) 
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Electrostatics 


This electric field is due entirely to the charge distribution, and so must be conservative, 
and we would expect that V x E = 0 as E is directed radially outward and so has no 
circulation. It follows that: 


V(r > a) 


V{r < a) 


Poa? 

(3 - (3)e 0 




poa 3 Ur') 2 3 /3(r') 1 "| r 

(3 - Wo 2-/3 + 3 . „ 


PQQ 3 

3(3-/3)(/3-2 )e 0 


(3 r 2 -^^" 3 


/3(/3 - 2 )r ^ . 


V(a)~ 


PQ 

3 £0 



r' dr ', 


(3 + 2/3 - /3 2 )p 0 « 2 
3(3 — /3)(/3 — 2)e 0 


Po 

3eo 


(r 


./\2 


(3 + 2/3 - /3 2 )p 0 a 2 _po_ 2 _ 2 
3(3 -/3)(/3 - 2)e 0 6e 0 ° 


(1.15) 

(1.16) 

(1.17) 

(1.18) 

(1.19) 

( 1 . 20 ) 
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Electrostatics 


1-3 The electric field is given by E(r) = Eq cos(z/zo) exp(—r/ro) r, where zo and ro are con¬ 
stants. Find the charge density. 

Solution 


In this problem the electric field is given in terms of z and r. We will need to write E in 
terms of either Cartesian or spherical coordinates, and then use Gauss’ law in differential 
form. Choosing spherical coordinates because E is in the radial direction, 


P( r) = £oV • E, 

= — 


r 2 cos 


r cos 6 
zo 


exp- 


( 1 . 21 ) 

( 1 . 22 ) 


e 0 Eo 


. r cos 0\ ( r\ 2 - ( r cos $ A cos 0 ( r , 

2r cos I - exp-— r sin - -exp-+ 

zo J \ ro J \ zo J zo \ ro 


r 2 cos 


r cos 6 
Zo 


eqEo ( r cos 6 

cos ( - I exp 

Zo 


exp 

r 

ro 


-) (- 

ro) \ r 0 

2 — tan 


r cos 6 \ r cos 6 r 


zo 


z o 


eo-^o f z\ ( r 

-cos — exp- 

r \z 0 J V r o 


^ i z \ z r 

2 — tan ( — )- 

Zo J z o r 0 


(1.23) 

(1.24) 

(1.25) 


1-4 If we had a point charge q at the origin we might choose the reference point to be some 
point at an arbitrary distance ro (usually infinity) from the origin. Then if we wish to 
find C(r, 0,0) it would be convenient to have the reference point at ro = (ro, #,</>). Al¬ 
though obtaining the potential in this case is trivial, and one would usually just write it 
down, obtain the potential by carrying out explicitly the line integral for an appropriately 
parameterised curve. 

Solution 
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We start by parameterising the path from ro to r: 

r'(A) (ro — A) r; dr' — — rfAr; (0 < A < ro — r). 

Then, 

V(r) = - f E(r')-dr', 

J 1*0 

rr(X=ro-r) 

= - / E(r'(A)) • dr 1 , 

Jr(. 


-I 

l 


r(A=0) 
r(A=ro— r) 


(A=0) 47T£ 0 (ro - A) 2 


r • (— dXr), 


q 


47re 0 (ro - A) 2 

1 r 0 -r 


d\, 


47T£o(ro - A) Jo 
Q 


q 


(1.26) 


47T£ 0 [?’o - (ro - r)] 47re 0 (ro - 0) ’ 

q _ q 

47T£or 47T£oro 

Hence, if we set ro = oo we get the usual potential for a point charge q at the origin 

q 


(1.27) 

(1.28) 

(1.29) 

(1.30) 

(1.31) 

(1.32) 

(1.33) 


V(r) = 


^TTSqT 


(1.34) 
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Electrostatics 


1-5 The electric field is given by E(r) = Eq cos(z/zo) exp(—r/ro) r, where zo and tq are con¬ 
stants. Check whether or not the electric field is conservative. If it is conservative find 
the potential, if it isn’t suggest how it may be possible to find the electrostatic part of the 
electric field (if present) and the corresponding electrostatic potential F(r). 

Solution 


First we need to test whether or not the field is purely electrostatic, i.e. whether or not 
it is conservative. If V x E = 0 then E is conservative. First write the field in spherical 
coordinates 


E(r, 0,0) = Eq cos(r cos 9/zo) exp(—r/ tq)t 


(1.35) 


and use 


V x A = 


r sin0 


5 dAg 

S ( " 6 A i)--^ 


^ 1 

r+ - 
r 


1 dAr _d_ 

sin# dcf> dr T ^ 


9 


1 

+ 

r 


^ rA>) 


dA r 

~W 


0. 


(1.36) 
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V x E = — 


ldE r 


r d6 


0, 


= -^exp(-r/r 0 )Jt cos(r cos 6/zq) <f>, 

Eq 

=-exp(—r/ro) [— sin(r cos9/zq)][— r sind/zo] <fi , 


Eq ^ 

-exp(—r/ro) sin(r cos 9/zq) sin 0 <fi. 

zo 


(1.37) 

(1.38) 

(1.39) 

(1.40) 


Since VxE/0 the electric field is not purely electrostatic. However, from Exercise 1—3 
we see that there is a non-zero charge density p(r, d, </>), and so there must be an electrostatic 
component of the electric field. This electrostatic field and potential could be computed 
from p using Coulomb’s law. 


1-6 How much work must be done to assemble: (a) a physical dipole made of charge +q and 
charge — q separated by distance d, (b) a physical quadrupole made up of four charges +g, 
—q, Eq and — q on successive corners of a square of side d, and (c) a physical quadrupole 
made up of four charges — g, +g, +q and — q equally spaced apart by distance d on a 
straight line (see diagram below). 


(a) 


• + 


d 


Solution 



(c) 


d. 


+ 


2 


+ 


3 


-m 

4 


The work done to bring together a group of N charges is 


1 N 

w = 


i=i 


(a) 


W = 


l EE + (+,) i^ 


ineod 


(1.41) 


(1.42) 
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(b) 


V 1 = 

qiV i = 

V 2 = 

42 Vi = 


: v 3 = 

q 1 

=i + ^ 

\ d d 

AtT£q ' 

43 Vi = 

4 2 

(V 2-2 y 

Aireod 

= U = 

q 1 

fl V2 

y.d d 

AtT£q ' 

44 V 4 = 

4 2 

(V 2-2 y 

Ansod 


Aireod 


2-2 


)• 


+ d 


47T£:od 


2 - V 2 , 


Hence, 


1 a" 5 

tU = - x 4 x ——- 

2 Aireod 


2-2 = 


2n£od 


(c) 


Vi 

Vi = 

Vi 

42 Vi = 


= Vi 

44 Vi 

= Vi 

43 Vi 


1 1 -1 

47T£o 2d 3<i 

4 


7 4 


6 47T£od’ 


7 ~ 2 

6 4tt£o^ 

4 


-1 1 -1 

47T£o v d d 2 d 

1 q 2 

2 47T£o^ 


1 4 


2 47T£orf’ 


Hence, 


W = 


4 


1 

_ ^ 

2 47T£o^ 


7 

6 


1 

2 


10 q 2 
3 Atteq d 


(1.43) 

(1.44) 

(1.45) 

(1.46) 

(1.47) 

(1.48) 

(1.49) 

(1.50) 

(1.51) 

(1.52) 


1-7 (a) Use Gauss’ law in integral form to find the electric field due to charge density p(r) = 
po exp(—r/ro), and (b) check that you obtain the original charge density by taking the 
divergence of the electric field you find. 
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Solution 


Gauss’ law 


E • dS = —Qenc, V • E = —. 
s £ 0 ^0 


( a ) 


47 ir 2 E r — — 




1 f r 

— / po exp(—r / /ro)47r(r / ) 2 dr , 

H) Jo 

r 0 (^-e ~ r '/ r °j (2rg + 2r 0 r' + (r') 2 ) 

-^2 2rg - r 0 e _r/Vo (2rg + 2r 0 r + r 2 ) 
eo?’ L 


^0 Jo 
Po 

r — o 


Jo 


(b) 


p = £qV • E 


(1.53) 


(1.54) 

(1.55) 

(1.56) 


(1.57) 
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1 d o 

= e ^l~ r 
r A dr 


Po 

r r .2 


~r o 


Po 
£o r 2 

-1 
^0 


= p 0 exp(-r/r 0 ). 


■ 2rg — roe r ^ r ° (2rg + 2ror + r 2 ) 
(2ro + 2r 0 r + r 2 ) - r 0 (2r 0 + 2r) 



(1.58) 

(1.59) 

(1.60) 


1-8 An isolated conducting sphere of radius a has net charge Q. Find how much work was 
done to charge the sphere using two different methods: (a) from the charge on the sphere 
and its potential, (b) by finding the energy stored in the electric field. 

Solution 


From Gauss’ law 


E r {r) 


0 (r < a) 

Q/Atts^t 2 (r > a) 


Q/Attsqci (r < a) 
(r > a ) 


(1.61) 


(a) The work done to bring together a group of N charges, or a continuous charge distri¬ 
bution p(r) is 


W = i / p(r)V(r).d 3 r. 

^ J all space 


(1.62) 


We can re-write this for a surface charge density 


\ f cr(r)V(r)dS , 

J s 

(1.63) 

IQ Q , 2 

2 4:7T£qCL 

(1.64) 

Q 2 

87TS0CL 

(1.65) 


(b) The energy density stored in an electric field is 


M r ) = y E ( r ) • E ( r )- 


( 1 . 66 ) 
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Since the sphere is conducting E = 0 for r < a, and so 


? / E 2 d?r, 

(1.67) 

z J all space 


2 L (w >) 4 " 2dr • 

(1.68) 

so q 2 A r ~ 2, 

n , 9 47T / r dr, 

2 (47re 0 ) 2 J a 

(1.69) 

Q 2 

(1.70) 
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Poisson's and Laplace's equations 


2 Poisson's and Laplace's equations 


2-1 Charge +g is located on the z axis a distance d/2 from a grounded plane conductor in 
the x-y plane. Find how much work was done to bring the charge to its current location 
using two different approaches: (a) the work done against the electrostatic force if the 
image charge were real and there was no grounded conductor, (b) the work done against 
the electrostatic force due to the induced surface charge 


( 2 . 1 ) 

where z is the height of the charge above the plane. 

Solution 


a(x,y,0) = 


27r (a; 2 + y 2 + z 2 ) 3 / 2 


(a) The force on charge +q at height +z due to image charge —q at height — z is 


F(z) = 

W(z = d/2) = 


1 


Q 


47T£o (2 z) 2 

(0,0, d/2) 


z. 


-L 


F • dr, 


(0,0,oo) 
2 


r 

Jdl 


z 2 dz, 


167T£o Jd/2 

1 t 

2 47reo d 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 



(b) From symmetry arguments, the force on charge +q at (0, 0, z) due to the real surface 
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charge density a(r) will be in the z direction, 


dF z (z ) 


= dF(z ) cos 6, 

_ (+g)[<7(r)d,S] 2 

47T£o R 2 R ’ 


(+<?) 

47 ts 0 R 2 2t r (i? 2 ) 3 / 2 i? 


9 * 4 2- dr, 


dF z (z) 

Fz(z) 

Fz(z) 

F(z) 


-gV 
4^0 i? 6 

"7 


r dr. 


z 2 r 


4:7T£o J o (r 2 + z 2 ) 3 


dr, 


~Q 


2 r 


4^0 

1 


4(r 2 + z 2 ) 2 
„2 


47T6:o ( 2z ) 2 


z. 


(2.7) 

(2.8) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


This is identical to the force on charge +q at height -\-z due to image charge — q, and so 
the work done will be identical to that calculated in part (a). 


2-2 Charge -j-q is brought near to two orthogonal grounded conducting planes, one correspond¬ 
ing to the x-z plane and the other to the y-z plane. The charge is located at (a, 6,0). 
Find the work done in bringing the charge from infinity to its current location (a) by using 
the method of images to find the potential at the location of the real charge, and (b) by 
considering the force on the charge as it is brought from infinity. 

Solution 


(a) At the location of charge -\-q the potential can be calculated as if it were due to the 
three image charges as in part (a) of the diagram below, 


V = 


47T£o 
N 


+ 


2a [(2a) 2 + (26) 2 ] 1 / 2 


W= 2 




2=1 


1 

2b 


r,; = 


167T£o 


1 1 1 

‘a + (a 2 + 6) 2 ] 1 / 2 ~b 


(2.14) 

(2.15) 
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Poisson's and Laplace's equations 


Note, it is only the real charge that enters into the sum above. 



x Q a x 


(b) We first calculate the force on charge +q at its final position due to the induced surface 
charge on the conductor as if it were due instead to the image charges as in part (a) of 
diagram above, 


F (a,b, z) 


F (a,b,z) 


ax 


+ 


by 


(+<?) ~gx +q 

4tt£o L(2a) 2 [(2a) 2 + (26) 2 ] V(a 2 + 6 2 )V 2 ' (a 2 + & 2 )V 2 


167T£o 


(a 2 + 6 2 ) 3 / 2 a 


-o x + 


(a 2 + 6 2 ) 3 / 2 


1 

62 


(2 b)\ 

( 2 . 16 ) 
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Poisson's and Laplace's equations 


Similarly, for the charge at some arbitrary position (x,y,z) the force is 


F(x,y,z) = 


X 


167T£o l\(x 2 + y 2 ) 3 / 2 x 2 


x + 


( x 2 + y 2 ) 3 / 2 y 2 


(2.17) 


The work done to move a charge from ri = (oo, oo, 0) to r 2 = (a, 6, 0) is W — — j^ F • dr, 


and because the electrostatic field is conservative, this is independent of the path taken. 
For convenience we split the path into two parts as in part (b) of the diagram above. Then 


W = - [ F(oo,y,0) • (-dy y) - [ F(x,6,0) • (-dxx), 

Jr! J r 2 


Q 


16 tt£o 
™2 


Ub 


OO / 1 \ PC 

-77>) dy + J a 




X 


-ft I dx 


(; X 2 + 6 2 ) 3 / 2 X 2 


Q 


167T£o 

J2 


Q 


1 

-T + 


+ 


1 


1 

+ - 


\fx 2 + 6 2 x\ a J ’ 


167T£o b (a 2 + 6) 1 / 2 


(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


which is the same as found in part (a). 


2-3 Show that the potential outside a long conducting cylinder of radius a in the presence of 
a long parallel line charge +A at distance d is identical to the potential of the line charge 
and a parallel image line charge —A at distance d\ from the cylinder’s axis towards the real 
line charge (see diagram below). [Hint: draw lines to point P from the two line charges. 
Use the cosine rule of triangles to write the two distances in terms of a, d and (/) and use 
the formula the for potential due to a line charge, and superposition, to write a formula 
for the potential at P. Finally require that V does not change if 0 changes.] 



Solution 
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Poisson's and Laplace's equations 



Using the cosine law: 


p\ — a 2 + d 2 — 2ad cos </>, p% — a z + df — 2 adi cos 


2 _ , ^2 


Adding the potentials at P of the real and image line charges, 


V(a,<j)) = 1 — [(+A) In pi + (-A)ln p 2 ], 


2tts 0 

1 




2 7T£ 0 \Pl 


1 1 n ( a 2 + $ — 2adi cos 6 

- In 


2tt£q 2 \ a 2 + d 2 — 2ad cos (f) 


For this to be constant on the cylinder’s surface, dV/d<fi = 0, i.e. 


d f a 2 + df — 2adi cos (j) 


2 adi sin 0 


d(f) \a 2 + d 2 — 2 ad cos (j) 

(a 2 + djf — 2adi cos (j))2ad sin 0 
(a 2 + d 2 — 2ad cos (j)) (a 2 + d 2 — 2ad cos </>) 2 

2a sin 0 [ di(a 2 + d 2 — 2adcos(j)) — (a 2 + d 2 — 2 adi cos <j))d\ 


= 0, 


= 0 . 


= 0 . 


(a 2 + d 2 — 2adcos(/)) 2 

(—d)d 2 + (d 2 H - a 2 )di -\- (—a 2 d) = 0. 


The solution of this quadratic equation is 


= a 2 /d or di — d. 

The physical solution is di = a 2 /d as di — d corresponds to the image line charge 
co-located with the real line charge +A — it is nevertheless a solution as V = 
cylinder’s surface, as well as everywhere else! 
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Poisson's and Laplace's equations 


2-4 Find the capacitance of a two-wire transmission line comprising two identical parallel cylin¬ 
drical conductors of radius a whose axes are separated by distance D (see diagram below). 
You may use the result for the potential due to a line charge near a single cylindrical 
conductor to find the potential difference by replacing the cylinders by equal but opposite 
image line charges, +A and —A (C m -1 ). The capacitance of two conductors with potential 
difference V and having charge +q on one and —q on the other is C — q/V 



The potential at A (and conductor 1 surface) due to the image line charges is, 
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Poisson's and Laplace's equations 


v A = - 


1 

2n£ 0 

[(+A) In pi + 

A 

In ( 


2n£ 0 

\ 

\P2 / 

A 

In 

f d — a A 

2tT£q 

\a-dij 


A _ ( d — a 
in 


2tt£o \a — a 2 /d 


——— In — 
2tt£q \a 


Similarly, the potential at B (and conductor 2 surface) due to image line-charge 


V B = + In f d ] . 


2tt£o \a y 

Hence the potential difference between the two conductors is 


Vba = (Vb-Va) = 

7T£ 0 \a' 


(2.31) 

(2.32) 

(2.33) 


—A is 


(2.34) 


(2.35) 


But d = D — di and di = a 2 /d, so 

.-. d= \{D + \Td 2 -4a 2 ) (2.36) 

is discarded because usually D a and 

d^> a. 

Hence the potential difference is 


d 2 — D d + a 2 = 0 (quadratic equation), 


since the other solution, d — ^(D — VD 2 — 4a 2 ) 


V BA = — In 

7T£0 




7T£o 


l^(D + VD* -4a 2 )^ 


(2.37) 


By definition the charge per unit length is A, and the capacitance per unit length (F m x ) 
is the charge per unit length divided by the potential difference, so that 


A 7V£q TTSq 

_ _ _ _ 

V BA In (^(D + V-D 2 - 4a 2 )) In (£>/“)’ 


(2.38) 


where the approximate result is valid for D a. 
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Essential Electromagnetism: Solutions 


Poisson's and Laplace's equations 


2-5 A region of space is bounded by three plane conductors as illustrated. Find the potential 
everywhere between the conductors. 


V(0,y)= 



Solution 


The potential must be finite at x — 0 and drop to zero as x —» oo, so we need the negative 
exponentials for the functions in x. At x = 0 the potential must be zero at y — 0 and so we 
need the sine functions for the functions in y. Furthermore V(0,b) — 0 requires k — rm/b 
so that the solution is 


v(x, y) = A n exp sin ■ (2.39) 

n =0 

The boundary conditions at x = 0 determine the coefficients A n in this Fourier sine series 


oo 



A 


n 


v Q . 



0 n even, 

4Vo/n7r n odd. 


(2.40) 

(2.41) 

(2.42) 

(2.43) 


2-6 Find the potential inside the rectangular region, 0<x< a, 0<y<b and 0 < z <c with 
V(x,y,c ) = Vq (x,y), and V—0 on the other 5 sides, where 


Vo(x,y) = V sin 



(2.44) 
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Poisson's and Laplace's equations 


Solution 


The potential for this case is of the form 


v ( X,y,z ) = EE A k i sinh ( 7 k i z) sin (a k x) sin (/?/ y ), 


k= 1 2=1 

kir Itt 9 2 ^2 

a k = Pi = — and 7 ^ = + p t . 


The potential at z = c may be written 


(2.45) 

(2.46) 


Vo (a;, y) = Vi sin (ana;) sin (foy) 


(2.47) 


and so the coefficients in the series are 


Am = 


4Vi 


n b 

sin {ol\x) sin (a k x ) sin (/? 3 y) sin (/?/ y) dx dy , (2.48) 

. 


a b sinh ( 7 ^ c) J 0 Jo 


and in this case there is only one non-zero coefficient 
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Poisson's and Laplace's equations 


^13 


^4l3 


v(x,y,z ) 


4Vi f a fb 

= —-—-- / sin 2 (a\x) dx / sin 2 (B^y) dy, 

a b sinh (713 c) j 0 V 1 ; J 0 KHdyj y ' 


4Vi 

x sin(2 a\x) 

a 

v 

y 

sin(2 0 3 y)~ 

a b sinh (713 c) 

2 4a 1 

A 

0 

2 

4/^3 


4Vl 

X 

sin(2 ttx i 

/«)] 

a 

X 

0 

y sin 

(67 xy/b)' 

a b sinh (713 c) 

2 

47 r/a 

- 

2 127 x/b 

4Vi 

a 

sin(27r) 

X 

~b 

sin(67r) 


a b sinh (713 c) 

2 

1 

2 ~ 

12vr/6 _ 



V ! 


sinh (713 c) ‘ 


= Vi 


sinh (7r/a) 2 + (37 t/6) 2 z 
sinh | ( 7 r/a) 2 + (37 t/6) 2 c 


/ 7VX\ 

sin ^—J sm 


3tt?/ 


(2.49) 

(2.50) 

(2.51) 

(2.52) 

(2.53) 

(2.54) 


Actually, we could have written down this answer straight away after recognising that 
Vo (x,y) was the product of one of the allowed functions of x having a = oq with one of 
the allowed functions of y having /? = ^3, from which we obtain immediately the solution 
in z with 7 = 713. 


2-7 The potential on a non-conducting sphere of radius a is given by 

V — Vo(3cos 2 9 + cos 9 — 1). (2.55) 

(a) Find the potential and electric field inside the sphere. 

(b) Find the potential and electric field outside the sphere. 

(c) Find the surface charge density on the sphere as a function of 9. 

Solution 

(a) Clearly we have spherical symmetry and no dependence on azimuthal coordinate </>. 
The general solution of Laplace’s equation with axial symmetry is 

00 

V(r, $,</>) = Z(A,r' + B< r (^ +1 )^ P^cos9). (2.56) 

e=o 
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Essential Electromagnetism: Solutions 


Poisson's and Laplace's equations 


Since the potential must be finite as r ^ 0 we must have = 0 for all t. Before applying 
the boundary condition it will simplify our working if we write it in terms of Legendre 
polynomials 

V(a,0,</>) = V o [Pi(cos0) + 2P 2 (cos0)]. (2.57) 


Then, applying the boundary condition, 

oo 

^2 At a e P e (cosd) = V o {Pi(cos0) + 2P 2 (cos0)], (2.58) 

£=0 

and by equating coefficients of Pg(cos9) we see that Ag = Vq/ a and A 2 = 2Vo/a 2 , giving 


V m (r,6,(l>) = 
V m (r,0,(p) = 
The electric field is 


-Pi (cos 9) + 2—^P 2 (cos9) 


- cos 9 H—^ (3 cos 2 9 — 1) 


Vo, 

V Q . 


(2.59) 

(2.60) 


. ,'dV^ ldV~ 1 9V-. 

E (r) — — r H -^v-0 4 :—> 

or r o9 rsin# o<p 


(2.61) 
V 0 0. (2.62) 


(b) Since the potential must tend to zero as r —> oo we must have Ag = 0 for all L Again, 
we write the boundary condition in terms of Legendre polynomials 


1 ^2 r, o . 


1 6r 

cos 6 H— 7T (3 cos 2 6 — 1) 

V 0 r + 

- sin 6 H— tt cos 6 sin 6 

a a 2 


a a z 


V(a, 9,4>) = K o [Pi(cos0) + 2P 2 (cos0)]. 


(2.63) 


Then, applying the boundary condition, 


Bga T +1 )p^(cos0) = Vo [Pi (cos#) + 2P2(cos0)], 
l=o 


(2.64) 
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Poisson's and Laplace's equations 


we 


see that B\ — a 2 Vo and £>2 = 2a 3 Vo, giving 


V out (r, 0 ,«/>) = 


v out M,</>) = 


r Pi(cos 0 ) + 2 ^-P 2 (cos 0 ) 


-2 cos 0 H—^ (3 cos 2 0 — 1 ) 


Vo, 

V 0 . 


(2.65) 

( 2 . 66 ) 


The electric field is 


\ ,'dV^ ldV~ 1 dV 

E (r) — — r H 0 H : v-wt 

1 or r o0 rsm 0 o<p 




(2.67) 


2 —^ cos 0 + 3 ^-( 3 cos 2 0 — 1 ) 


V 0 f + 


v a? 


-it sin 0 + 6—r cos 0 sin 0 


V o 0. 


( 2 . 68 ) 


(c) We use Gauss’ law in integral form for a small section of the sphere of area 8S located 
at (a, 0 , <j>) inside an infinitesimally thin gaussian pill box having an upper surface of area 
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Essential Electromagnetism: Solutions 


Poisson's and Laplace's equations 


5S just outside the sphere and a lower surface just inside the sphere. For the upper surface 
of the pill box the normal unit vector outwards from the pill box is n = r , whereas for the 
lower surface of the pill box the normal unit vector outwards from the pill box is n = —r. 
Applying Gauss law in integral form 

E out (a,6,(t>) ■ (5Sr) +E in (a,e,(f>) ■ (~5Sr) = a{6)5S/£ 0 . (2.69) 

Hence, 

cr(0) =e 0 [E^\a,d^) (2-70) 

= £o— ([2cos# + 3(3cos 2 6 — 1)] + [cos# + 2(3cos 2 # — 1)]) , (2.71) 

= £ 0 ^ (3 cos # + 15 cos 2 # - 5) . (2.72) 


2-8 Consider a point charge on the z-axis at z = r f . Find V(r, #,</>) in terms of Legendre 
polynomials for r > r '. 

Solution 

There is no dependence of V on 0, so 

00 

V(r,e,<t>) = J2( A nr n + B n r-^)P n (cos6) 

71=0 


(2.73) 


The diagram shows the geometry for this problem. 



The boundary condition for this problem will be the potential along the z axis for z > r'. 
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Poisson's and Laplace's equations 


which we can obtain from Coulomb’s law 


V(r, 0,(f>) = j^—(r - r') \ 

4:7T£o 


Q -l 
r 


4:7V£q 


r 


°° f T f \ n 

- -r _1 ( — ) (binomial series). 

47r£ ° P 0 w 


(2.74) 

(2.75) 

(2.76) 


On the z-axis 6 — 0 so P n (cos6) = P n ( 1) = 1, and the general solution for the potential is 


00 

V(r, 0 A) = Y j (A n r n + B n r~^ n+1 ^ 


(2.77) 


n =0 


Hence, 


J2(A n r n + B^- 1 ) 


n —0 


00 / / \ n 


(2.78) 


giving 


— 0 ) Pn — 


( r') n q 
47nm ’ 


V(r > r',9,(j)) = - y^(r') n r ( n + x )p f cos m 

47r£o 


(2.79) 

(2.80) 


77=0 


2-9 The potential on the surface of a sphere is 

V (a, 0, 0) = V\ sin 6 sin </> + V 2 sin 0 cos 0 sin </>. (2.81) 

Find the potential inside the sphere. 

Solution 

The solution is of the form 

00 t 

V(r,0,ct>) = J2 J2 + ( 2 - 82 ) 

1=0 m =—£ 
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Poisson's and Laplace's equations 


The requirement that T^(0, 0, (j>) is finite gives Bi = 0 for all £. The boundary condition at 
r — a can be re-written in terms of spherical harmonics as follows 

I O '!r f O7T 

y(a,0,^) = Uy y [U,i(M) - U,-i(M)] + v 2 ^j— [r 2 ,i(M) ~ U,-i(M)] . 

(2.83) 

(2.84) 

(2.85) 

giving 


Hence, comparing coefficients we find 



V(r,e,<l>) 


2 

V\ - sin 9 cos (f) + V 2 ^ sin 9 cos 9 cos (j). 
a a z 


[Y 2 , 1 (e,<f>)-Y 2 ,-i(e, </>)}, 
( 2 . 86 ) 
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Multipole expansion for localised charge distribution 


3 Multipole expansion for localised charge distribution 


3-1 On the surface of a non-conducting sphere of radius a is surface charge density cr(a, 0, (j)) 
ao cos 3 9. Find the dipole moment of the sphere. 

Solution 



For a surface charge density cr(r) the dipole moment is 

p = J a( r) r dS. 

In this example, the surface charge density depends only on 9 and so p = p z z where 


(3.1) 


Pz 


= J cr( r) z dS 

r 7 T 

= / cf{9) z 2na 2 sin 9 d9, 

J o 

= J (cro cos 3 9) (a cos 9) 2na 2 d(cos0), 


= 2tvcl 3 <jo 

4 3 

p z = -7 TCL CT 0 . 

5 


P = -7 TCL 

5 


cos 5 9 


5 nl 1 


-1 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


32 


Download free eBooks at bookboon.com 















Essential Electromagnetism: Solutions 
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3-2 The quadrupole potential is 


Vuad(r) - 47r£or 3 / p{v')(r'f 2 [3(r • r ') 2 l] d\'. 

(3.8) 

Show that it can be written as 


i i 3 3 

tquadM - 4 r 5 2 EE^ r ’^ 

1=1 j = 1 

(3.9) 

where the quadrupole moment tensor is 


Qij = J p( r ') [3 r'r- - Sij (r') 2 ] dV. 

(3.10) 

[This exercise is easy using index notation.] 


Solution 

The quadrupole potential is 


Kuad(r) = 47r£or 3 J p(r')(r') 2 2 [3(r • r ') 2 l] dV, 

(3.11) 

= 4t r£ 0 r 5 / p{V ' )r2 (/)2 2 ^ ‘ r ,)2 ^ dV ’ 

(3.12) 

= 4t r£ 0 r 5 j ^ } 2 t 3(r ‘ r )2 ^ (r )2 ] ^ ’ 

(3.13) 

= 47r£or 5 / ^( r ) 2 [ 3 rir * r Vi rir * ( r ^ ’ 

(3.14) 

= 4t r£ 0 r 5 / P(r )r ' rj 2 ^ (r )2 ] ^ ’ 

(3.15) 

Viuad(r) - 47r£or 5 2 ^ r?r > 

(3.16) 


Index notation and the Einstein summation convention has been used above, but writing 
the summation explicitly we have 

1 i 3 3 

Vuad(r) = 47r£or 5 2 X X Qv rir i ( 3 - 17 ) 
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Multipole expansion for localised charge distribution 


where 


Qij 



(3.18) 


3-3 A physical quadrupole is made up of four charges lined up along the 2 : axis: -qo at (0,0,—2a), 
+go at (0,0,—a), +go at (0,0,a) and -qo at (0,0,2a). (a) Obtain the quadrupole moment, 
(b) Find the potential at r = (6, 6, 0) for b >> a. 

Solution 


charge #12 3 4 

-q +q +q -q 

-O- • - 1 - # -• 

-2a -a 0 a 2a 

(a) The quadrupole moment tensor is 


Qij — 



Sij (r') 2 ] d 3 r'. 


z 


(3.19) 
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Multipole expansion for localised charge distribution 


For N point charges this becomes 


N 

Qij = ^ ^ Qk 

k =1 


3 - 5ij (r^Y 


(3.20) 


where r^ is postion of charge q^. 

Since r{^ = rj^ = 0 for all four charges as they are on the z axis, the quadrupole moment 
tensor is diagonal with only Q n, Q 22 and Q 33 being non-zero, 


4 

Qll =^2vk 

k =1 


3x0 x 0-(5 n (r [fc1 ) 2 , 


(3.21) 


— ( -< ?) [-(- 2 a) 2 ] + (+?) [ _ ( _a ) 2 ] + (+?) [ _ ( a ) 2 ] + (“<?)[“( 2 a) 2 ] , 


Qn = 6 qa 2 . 


(3.22) 

(3.23) 


Similarly, 


Q 22 = 6 qa 2 . 

Q33 = ^2qk 3 x rf ] x rf ] - <5 33 (r [/c] ) 2 


k =1 


= (- 9 ) [3(-2a ) 2 - (2a) 2 ] + (+q) [3(-a ) 2 - (a) 2 ] 

+ (+<?) [3(a ) 2 — («) 2 ] + (— q ) [3(2a ) 2 — ( 2 a) 2 ] , 

Q 33 = — 12 ga 2 . 

6 qa 2 0 0 


Qij — 


0 6 qa 2 0 

0 0 — 12 ga 2 


(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 


(b) The potential for r j$> a can be approximated by the quadrupole potential 


1 1 3 3 

^quad(r) = 4 r 5 | EE Qw- 

1=1 j = 1 


(3.29) 


Only Qn , Q 22 and Q 33 are non-zero, and for r = (6,6,0) the distance from the origin is 
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Multipole expansion for localised charge distribution 


r = y/2.b, so that 


VquadO, &, 0) 


1 1 

47 T£o(v / 2 fr ) 5 2 


[ 6 qa 2 x b x b + 6 qa 2 x b x b 


3 qa 2 

87t\/2^o^ 3 


12ga 2 x 0 x 0] , 


(3.30) 

(3.31) 


3-4 Charge — q is located at the origin and charge is located at 

(x,y,z) = (a sin 0 o cos 0 o, a sin 0 o sin 0 o, a cos 0 o)- 

(a) Find the the non-zero moments of the multipole expansion of the potential in Cartesian 
coordinates, i.e. g, p, Qij (if non-zero), and use these moments in the multipole expansion 
in Cartesian coordinates to find the potential at (x, y, z) = (r sin 0 cos 0 , r sin 0 sin 0 , r cos 0) 
where r a. 

(b) Find the non-zero moments of the multipole expansion of the potential in spherical 
coordinates, i.e. 

Qirn = f Yl m {e'A')T d p{v')d\', (3.32) 

and use these moments in the multipole expansion in spherical coordinates to find the po¬ 
tential at (r, 0, 0 ) where r a. Compare the result with that from part (a). 

Solution 


The net charge (monopole moment) is zero. There are two equal but opposite charges and 
so we have an electric dipole moment, and no higher moments. 

(a) In Cartesian coordinates 

p = (—q) (0,0,0) + (+g) (a sin 0o cos 0o, a sin 0 q sin 0o, a cos 0o), (3.33) 

— Po x (sin0o cos0o, sin0o sin0o, cos0o), (3.34) 

where po = qa- The potential at (x,y,z) = (rsin0 cos0,rsin0 sin0,rcos0) where r a 
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is 


V(r) = 


4:7T£o r 

Po 

4:7T£Qr 

Po 


^p-r, (3.35) 

2 (sin cos 0o, sin 0o sin 0o, cos ) • (sin 0 cos 0, sin 0 sin 0, cos 0), (3.36) 


47T£o r 2 
Po 

4:7T£o r 
Po 

4:7V£q r 2 


(sin 0o cos 0o sin 0 cos 0 + sin 0 q sin 0o sin 0 sin 0 + cos 0o cos 0), (3.37) 
2 [sin 0o sin 0 (cos 0o cos 0 + sin 0o sin 0) + cos 0o cos 0] , (3.38) 


[sin 0o sin 0 cos(0o — 0) + cos 0o cos 0] 


(3.39) 


(a) In spherical coordinates the multipole moments are given by 


I Y e * m (9',<j>')r re p( r')dV. 


(3.40) 


Writing the charge density using Dirac delta functions in spherical coordinates we have 


> Apply now 
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Multipole expansion for localised charge distribution 


P( r ) = ~qS(r) +q 


6{r - a)5(6 - 0o)5(<f> - fa) 
r 2 sin 9 


(3-41) 


There will only be dipole (t — 1) multipole moments 


r27T /*7T poo 


/*7T pC 

J 0 J 0 



. aid, qS(r - a)5(9 - 6 0 )5(4> - 0 0 ) 2 . 

qi-i — I I I \ I — sin 9e 9 r -~--- -r smOdrdOdd, 

87r r z sinO 


9i-i = aq\ — sin 

o7T 


(3.42) 


r2n rn r oo / ^ 

Jo Jo Jo V 4vr 


910 


cos r ^(r-«W»-».W-7.) r 2 8in edr ie dlp 
r z sin 0 


9io = aq\ — cos/9,). 
4-7T 


(3.43) 


9ii 


p 2n p 7T /‘OO 

/o Jo Jo 



sin 9e -i* r ain Wr ^ 

87 t H sin 0 


9ii = -09l/ ysin0 o e 

57T 


(3.44) 


The potential for r a will then be the same as in part (a), 


oo £ 


v(r,e,<j)) = (€+1)y i?,m(M), 


^=0 m——i 


(3.45) 


aq 

3eor 2 



sin 


r 3 _ i( k I 3 . / 3 

— sin 0e 9 + W — cos 0q\ — cos 0 
87T V 47T V 47T 


+ \ —sin 9 0 e-^°\ — sin de** 
\ 8n \ 8n 


(3.46) 


aq 

2>£ 0 r 2 


aq 

Attsqt 2 


-f-sin 0o sin 0e^° ^ + -^-cos0ocos0 
87T 47T 


H-sin 0o sin 0e ^ 

87T 


[sin 0o sin 0 cos(</>o — </>) + cos 0o cos 0] . 


(3.47) 

(3.48) 
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Macroscopic and microspcopic dielectric theory 


4 Macroscopic and microscopic dielectric theory 


4-1 A dielectric sphere (dielectric constant K) of radius a is placed in an initially uniform 
electric field Eo. (a) What are the boundary conditions on V , E and D for this problem, 
(b) Find the potential everywhere, (c) Find E, D and P everywhere, (d) Find the dipole 
moment of the sphere and the surface polarisation charge density. 

Solution 


(a) The boundary conditions at the surface of the sphere are that E ||, D± and V are 
continuous across the boundary. In addition the electric field very far from the sphere 
must equal the initial field. Defining this to be in the ^-direction, 

E(r » a, 9 , </>) = E 0 z, (4.1) 

.'. V(r a, 0, (f)) — — Eqz — — EqtcosO = — EbrPi (cos 0). (4.2) 

Since the potential has not been specified anywhere, we are free for convenience to set 

v(o,M) = o. 

(b) This is a problem with spherical symmetry but with no dependence on <j). Hence, we 
can write down the form of the potential 

OO 

V(r,0,4>) = ^ (A e r e + £<?r _(m) ) P e (cosO). (4.3) 

e=o 

So that Vi n is finite inside the sphere (containing r — 0), we must have all B™ = 0. 
Similarly, so that V ou t is finite outside the sphere we must have all A^ ut = 0, except as 
needed to give V{r^> a, 0, </>) = — EorPi(cos0), i.e. A° ut = — Eq. Hence 


Vin(r, 0, </>) = Af r e Pt( cos 0), (4.4) 

e=o 

oo 

Vut (r,0,<t>) = -E o r 1 P 1 (cos0) + ^ r~^P e (cos0). (4.5) 

z=o 


Applying the boundary condition on V at r = a, and remembering that we set V (0, 0, (j>) = 
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0, 


oo oo 

J2 A f a£p e( cos °) = ~E 0 a 1 P 1 (cos9) + J2 B r ta ~ (£+1)p i( cos#), (4.6) 

e=o e=o 

Afa = -E 0 a + B° ut a~ 2 . (4.7) 

with all other coefficients zero. Hence, we can now write the form of the solution as 

v in (r, 9, (j)) = Afr cos 9, (4.8) 

V ou t{r,9,(f)) = — Eqv cos 9 + B° ut r~ 2 cos 9. (4.9) 


Next we apply the boundary conditions on E and D at r = a, 


E|l = 

A'™ sin 9 = 


1 dV_ 

r d9 


0 , 


-E() sin 9 + B° ut a 3 sin 9 


(4.10) 

(4.11) 
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_ 8V 


—eA l i cos 9 = £qEo cos 9 + 2s q B° ut a 3 cos 9. 


Equations 4.7, 4.11 and 4.13 are three equations in two unknowns, but we only 
equations. Solving Eqs. 4.11 and 4.13 gives 


= E oa 3 7- 0 £ °\ , 4 n = -Eo 3e ° 


(s + 2^o) 


(£ + 2^o) 


Hence, 

MnM,0) = 

Vout {r,9,4>) = 


^ 3g() /) 77 3g() 

E/q / _ i 0 x r cos 9 —Eq f _ , o N z, 


(£ + 2^q) 

—Eor cos# + Eoa' 


3 (g — gp) ^_ 2 
(g + 2go) 


(g + 2go) 
r ~ 2 cos#. 


(c) The electric field is E = -W and the displacement field is D = gE 


dV 1 8V 
E(r,0,0) = -—r-- —0 


1 -7 

~WT V- 


Ei n — 


dr r 89 r sin0 8(f) 
3go 


(g + 2go) 


Dp, Din — ^Dj n . 


Eout = Eo + /T 0 £Q \ a 3 Eo 

(g + 2g 0 ) 


2 cos 0 r + sin 9 6 


T* j D out — ggEout- 


The polarisation field inside the dielectric is obtained from P = D — gpE giving 




(d) Since the polarisation field is uniform the dipole moment is 


4 o 4 ,3(s — ecO^o,-, 

P = o™ P = o™ 7 , n N E 0 . 

3 3 (g + 2go) 
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The surface polarisation charge density on the sphere is 


cr p (a, 0, (j>) = P • n = 


3(g ~ gp)gQ 

(s + 2^o) 


Eq cos 9. 


(4.22) 


4-2 An electret, i.e. a piece of material with a permanent electric polarisation, is in the shape 
of a sphere of radius a and has P(r) = Pq. (a) Find the surface polarisation charge density 
and the dipole moment of the sphere, (b) find V, E, and D everywhere, and (c) sketch the 
field lines of E and D. 

Solution 


(a) We are free to choose the sphere to be polarised in the z direction. Then 


^poi(<T$ 5 0) = P n = (P 0 z)-r = P 0 cos 9 = PqPi(cos9). (4.23) 

As the polarisation is uniform, we can obtain the dipole moment directly from P and the 
sphere’s volume 


P = 


7ra 3 Pn 


z. 


(4.24) 


(b) We first need to write down the form of the solution for the potential (inside and outside 
the sphere), and then apply the boundary conditions to fix the coefficients in the series for 
V. 

The form of the solution for the potential is 


V(r, 0,0) = + B^r (^+ 1 ))p^(cos#). 

t=o 


(4.25) 


Examining the angular dependence of cr po i we realise that the solution will only involve 
terms with £ < 1, then 


E m (r, 0 , (/)) = Ao + Air cos 9, (4.26) 

E out (r,0,0) = Bqt -1 + B\r~ 2 cos 9. (4.27) 
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The potential must be continuous ar r = a so that 

A 0 = B 0 a~ 1 , A 1 = B 1 a~ 3 . (4.28) 

Gauss’ law can be used to provide a boundary condition on E, and for this we will need 
the normal (in this case radial) components of the electric field, E r = — dV/dr , 

i^ n (r, 0, (/)) = —Ai cos 0, (4.29) 

£'° ut (r, 9, <f>) = B 0 r~ 2 + 2B ir ~ 3 cos 0. (4.30) 

Gauss’ law applied to the a pillbox spanning r = a at (a, 0, 0) is then, 

£'° ut (a, 0, (j)) - E^(a, 9,4>) = <r po i(a, 9, <f>)/e 0 , (4.31) 

Boa~ 2 + 2I?ia -3 cos 9 + A\ cos 9 = Pq cos 9/e q. (4.32) 


From this we see that Bq = 0, and then from Eq. 4.28 that Aq = 0, and that 
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Ai = Po/3^o, B\ — a 3 Po/3^o. 


(4.33) 


Hence, the potential is 


V in (r,6,<P) 

V ou \r,0,cf>) 


Po a Po 
-—r cost) = - — z, 
3eo 3eo 

a3p o -2 a 
- r cosy. 

3e 0 


The electric field is 


E(r,M) 
E in (r,M) 
E out (r, 0, <f>) 


8V ^ _ 1 8V - _ 
dr r 80 

_P^~ = P^ 
3eo 3eo 


1 dV 4> 

r sin# d(f> 

(— cos0r + sin 06 



2cosy? + siny^ , 


(4.34) 

(4.35) 


(4.36) 

(4.37) 

(4.38) 


The electric potential and field outside the sphere is identical to that of a dipole 


1 p-r _ 1 p cos 6 

dip 4:7tsq r 2 4:7 tsq r 2 5 

2cos0r + sin 0 0^ . 

and is consistent with p calculated earlier from Pg and the volume. 
Finally, the displacement field is given by D = £qE + P, 

D in M,0)Pp o , 

D out (r, y, <f>) = - ^ (> r~ 3 ^2cosyr + siny0^ . 


Edip(r) = 


P 


4:7T£q r 3 


(4.39) 

(4.40) 


(4.41) 

(4.42) 


(c) To sketch the electric field lines we notice that inside the sphere E is constant and in 
the —z direction, and that at the poles |P m | = |P out |/2, and that outside the sphere it has 
a dipole field. Also, electric field lines start on positive charge (either free or polarisation 
charge) and end on negative charge (free or polarisation charge). 

To sketch the displacement field we notice that inside the sphere D is constant and in the 
+z direction, at the poles that |P m | = |P out |, and that outside the sphere it has a dipole 
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field. Also, since there is no free charge present the field lines of D must form closed loops. 



4-3 The space between two concentric conducting cylinders of radius a and b > a and length 
L b is filled with a dielectric with permittivity e. The inner and outer conductors are 
held at potentials V a and V 5 , respectively. Find: (a) E, D and P everywhere; (b) the po¬ 
larisation surface and volume charge density everywhere, and the net polarisation charge; 
(c) the free charge on the inner and outer conductors, and the capacitance. 


Solution 



(a) We start by solving Laplace’s equation in cylindrical coordinates with no dependence 
on (j> and z 


1 d_ ( dV\ 

p dp \ dp ) 


(4.43) 


Integrating, we get 


dV 


P 


dp 


= A, 


J dV — A 


dp 

P 


V = A\np + B. 


(4.44) 
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A and B are integration constants to be determined from the boundary conditions at p — a 
and p — 6, 


V a = A In a + B, = A In b + B. 


Solving for A and B , 


„ 04 -v a ) D T/ , (14 - V b ) 

A = \ /T , / , B = V a + ; ~—;—~~ In CL. 


In (6/a) 

^(p) = T4 + (14 -14) 


ln(6/a) 


In(p/a) 

ln(6/n) 


(4.45) 


(4.46) 

(4.47) 


The electric field will be present only between the inner and outer conductors 


E(p) = - 


dV, 
dp 1 


(Vb ~ Va) V 

In (b/a) p l 


(4.48) 


Since the dielectric is linear the displacement and polarisation fields, again only between 
the inner and outer conductors, are 



Plug into The Power of Knowledge Engineering. 
Visit us at www.skf.com/knowledge 


Brain power 


By 2020, wind could provide one-tenth of our planet’s 
electricity needs. Already today, SKF’s innovative know¬ 
how is crucial to running a large proportion of the 
world’s wind turbines. 

Up to 25 1 of the generating costs relate to mainte¬ 
nance. These can be reduced dramatically thanks to our 
^sterns for on-line condition monitoring and automatic 
lulhkation. We help make it more economical to create 
cleanSkdneaper energy out of thin air. 

By snlkcjour experience, expertise, and creativity, 
industries can boost performance beyond expectations. 

Therefore we need the best employees who can 
ineet this challenge! 


The Power of Knowledge Engineering 


46 


^3 


Click on the ad to read more 


Download free eBooks at bookboon.com 






Essential Electromagnetism: Solutions 


Macroscopic and microspcopic dielectric theory 


n _ p = AVa-n) i- 

In (b/a) p 

r» 771 , ^(Va-Vb) 1 ^ 

P = d_£ oE = (e - go) - p. 

m (b/a) p 


(4.49) 

(4.50) 


(b) The volume polarisation charge density is 


Ppoi — v • p — ^ dp(pPp) ~ 0- 


(4.51) 


The surface polarisation charge density is cr po i = P • n, 

,7 , ,(V a -V b )l^^ , (V a -V b )l 

O’pol(a) = (g-£o) ^JT7T\ ZP P = ( £ ~ £ 0 ) 


ln(6/a) a 


m b (Va - V b ) 1^ ( ^ ^ 

^0,(6) = (e-co)-^-- b P (-P) = -(*-*0) 


In (b/a) a 

(K ~ Vb) 1 


In (b/a) b 


(4.52) 

(4.53) 


The net polarisation charge is 


Qpol = L X [27ra X cr po i(a) + 2irb X cr po i(6)] , 

Tt ,(V a -V b ) ( 2na 2irb\ 

= L{e - eo) -HbJaP{—-—) = °' 


(4.54) 

(4.55) 


(c) We obtain the free charge present on the conductors using Gauss’ law from which 
a/ = D • n, 


°f( a ) 

°f( b ) 


(' V a -V h )l _ (V a -v b )l , , iy a -v b ) /A ^ 

e , /T , - pp = e ——-’ 9f\ a ) = 27rLg 1 fuj - , (4.56) 


6: 


ln(6/a) a 

(V a ~ V b ) 1 


ln(b/a) a 


ln(6/a) b 


-- , ^ 04 - 14 ) 1 m 0 T 

rP’(-p) = ~ g i~Tru7Z\ T~? ?/(&) = -2ttL£: 


ln(6/a) 

(K - H) 


ln(6/a) 6’ 


ln(6/a) 

(4.57) 


Hence, the capacitance is C = qf/{V a — V&), 


C = 


2tiLs 
In (5/a)’ 


(4.58) 
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4-4 A spherical capacitor is filled with two different dielectrics with permittivities S\ and e 2 as 
shown in the diagram. The capacitor is charged such that charge is on the inner con¬ 
ductor. Find: (a) D, E and P everywhere; (b) the polarisation surface and volume charge 
density everywhere; (c) the net polarisation charge; (d) the potential difference between 
the inner and outer conductor, and the capacitance of the capacitor. 



Solution 


Because of the spherical symmetry we can use Gauss’ law in integral form to find the 
displacement field between the conductors 


47 xr 2 D r = +q, D(a<r <c, 9, 6) = -— 

47 xr z 


The electric field is E = D/e, 


E(r,M) = 


(a < r < b) 
( b<r<c) 


The polarisation field is P = D — £qE, 


P (r,M) = 


Q 

47rr 2 
Q 

47rr 2 


(a < r < b) 
(b < r < c) 


(4.59) 


(4.60) 


(4.61) 


(b) The volume polarisation charge density is 

Ppo\ -V P = - ^ ^( r2p r) = 0. (4.62) 

The surface polarisation charge density is a po \ = P • n, and each of the two dielectrics will 
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have surface polarisation charge at its inner and outer radii. Hence, 






pol 




"Sw = 


pol 


47ra 2 

q 

4:7Tb 2 

q 

47rb 2 

q 

4ttc 2 








^2 . 


1-^1 

47ra 2 V ei 1 ' 


1-^ 


47t6 2 

q 


Si 


, 1--1 

4tt6 2 V S2 1, 


1 - — 


4ttc 2 


s 2 


(c) The net polarisation charge is 


g po i = 4 tt a 2 a^{a) + 4t rb 2 [(7^(6) + <tS(&)1 + ^c 2 a { ^ x (c)i 


= q 


£i 


£l 




£2 


£2 


= 0. 


(4.63) 

(4.64) 

(4.65) 

(4.66) 

(4.67) 

(4.68) 
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(d) The potential difference is V a — V c = — J^ E • dr, 


V b -V c 

V a ~V b 

V a ~V c 


47T E2T 2 
Q 


dr = 


--f 

= -rj-s-j* = [ 

Jb 4vreir 2 J b 


Q 


c A.TlE2r 2 
1 q 


z dr = 


Q 


4:7T£ir 2 


dr = 


'1 _ 1 
47r^2 V b c 

1 1 


47T6:i V a b 


47T 


^ /I _ 1\ ^ /I _ 1 

la 6/ £2 \b c 


Hence, the capacitance is 


C = 47t 


_1_ /I _ 1\ _1_ /I _ 1 

8 \ \a bJ 82 \b c 


i -l 


(4.69) 

(4.70) 

(4.71) 


(4.72) 


4-5 A uniform slab of material with permittivity e\ is suspended parallel to the rry-plane, and 
has its lower surface at z = 0 and its upper surface at z — d. Outside the slab there is 
a uniform electric field Eo = £7>(sin0ox — cos 9q z). (a) Find formulae for E, D and P 
in the dielectric, the angle between E in the dielectric and the normal to the surface, and 
the surface polarisation charge density at z = 0 and z — d. (b) Find numerical values for 
the case of Eq — 1000 V m -1 , 9q — 45° and e r = 2, and include a sketch showing field 
directions. 

Solution 

A field line will bend as in the diagram. 



(a) The component of E parallel to the boundary is unchanged, and since there is no free 
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charge the component of D normal to the boundary is unchanged 

Eq sin 9q — Ei sin 9\ , £qEq cos 0o = e\E\ cos 6\ . (4.73) 

The electric field is 


Ei Ei sin 0ix — Ei cos 0iz. 

Ei = Eq sin 0ox — —Eq cos 9qz. 

£1 

Hence, from Eqs. 4.75 the magnitude and direction 


Ei 


sin 2 0o + 



11/2 


cos 2 00 


£0, 


01 — arctan 


( — tan 0o 

V^o 


(4.74) 

(4.75) 


(4.76) 


The displacement and polarisation fileds are 

Di = £iEi = £iEosin0ox — £qEqcos9qz^ (4.77) 

Pi = (Di-e 0 Ei) = (ei - g 0 )P 0 sin d 0 £ - ^ ~ £ °' >g ° P 0 cos 6> 0 z. (4.78) 

£l 

The surface polarisation charge is 

Cpoi(®,y,0) = +P\ cos 6 \, o-poi (x,y,d) = -Pi cosi^. (4.79) 

(b) Substituting for the case of Eq = 1000 V m -1 , 6q = 45° and e r = 2, i.e. e\ = 2eo we 
find 


01 = 63.4°, Pi = 791 Y m -1 , Pi = 1.40 x 10“ 8 C m" 2 , (4.80) 

p l = 7.00 x 10“ 9 C nr" 2 , (4.81) 

a po i(x,y, z=0) = +3.13 x 10 -9 C m -2 , (4.82) 

<Tpoi(x, y,z = d) = -3.13 x 10“ 9 C m" 2 . (4.83) 


4-6 Derive the force on an electric dipole in a non-uniform electric field. 
Solution 
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Consider a physical dipole consisting of charge +q located at r pos and charge — q located 
at r neg . Its dipole moment is p = qd where d = (r pos — r neg ). It follows that the force is 


F — (■P < ?)E'( r pos) "i" ( (?)E( r neg)> 

(4.84) 

= gAE 

(4.85) 

= <?(x A E x + y A E y + z A E z ), 

(4.86) 

= q [x(d • VE X ) + y (d • VE y ) + z (d • VE Z )\, 

(4.87) 

= q [(d • VE x x) + (d • VE y y) + (d • VE Z z)\, 

(4.88) 

= ?(d-VE), 

(4.89) 

= (qd ■ V)E, 

(4.90) 

F = (p • V)E. 

(4.91) 
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4-7 The relative permittivities of Nitrogen, Argon and Hydrogen in gas (at 20°C) and liquid 
phases are given below. 


Element 

N 2 gas 

Ar gas 

H 2 gas 

N 2 liquid 

Ar liquid 

H 2 liquid 

E r 

1.000546 

1.000517 

1.000272 

1.45 

1.53 

1.22 


http: / / www.kayelaby.npl.co.uk/general_physics / 2_2 / 2_2_1 .html 


Use the Clausius-Mossotti formula to find the electronic polarisability, and compare the 
results for the same elements in the liquid and gas phases. [You will need to look up any 
constants and the atomic weights and densities required.] 

Solution 

We need to use the Clausius-Mossotti formula 


3^o ( e r - 1 \ 
" pol_ N U + 2 J 


(4.92) 


where N = p/(Au ), p is the density, A is the mean atomic mass, u = 1.66 x 10 27 kg and 
£ 0 = 8.85 x 10“ 12 F m -1 . 

Values of density in the liquid and gas phases, and the mean molecular weight have been 
looked up in tables of physical/chemical constants and have been added to the table, and 
the polarisability calculated using Eq. 4.92. 


Element 

N 2 gas 

Ar gas 

H 2 gas 

N 2 liquid 

Ar liquid 

H 2 liquid 

£ r 

1.000546 

1.000517 

1.000272 

1.45 

1.53 

1.22 

p (kg m -3 ) 

1.25 

1.78 

0.089 

800. 

1393. 

67.8 

A 

14. 

39.95 

1 . 

14. 

39.95 

1 . 

a pol (C m 2 V -1 ) 

8.98 xl0“ 41 

1.70 xl0“ 4 ° 

4.49 x 10 -41 

1.01 xlO“ 40 

1.89 x 10 -40 

4.44 x 10 -41 


The agreement in <a po i for the same element between liquid and gas phases is quite good, 
the difference being at most 10%. 
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5 Magnetic field and vector potential 


5-1 Two parallel wires are separated by distance a and carry currents I\ and I 2 in the same 
direction. Find the force per unit length of wire. Include a diagram showing the direction 
of the force. If I\ — I 2 = 1 A and a — 1 m, what is the magnitude of the force per unit 
length? 

Solution 



The force on a circuit in a magnetic field is given by Ampere’s force law 


Fmag = j)(Idr X B). 


(5.1) 


For parallel currents as in the diagram, the force on length L of wire 2 is 


F — I 2 LB 1 = I 2 L 


2na 


(5.2) 


and is directed towards wire 1 as shown. The force per unit length is 

F_ _ Ushh 

L 2iva 


If the wires are 1 m apart and each carry 1 A, 


F 4i r x 10 


-7 


27r 


= 2 x 10“ 7 N m" 1 . 


(5.3) 


(5.4) 


The amp is defined as the current flowing in two parallel wires 1 m apart such that the 
force beween them per unit length is 2 x 10 -7 N m -1 . 
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5-2 Using the equation for the vector potential in terms of the current density find the vector 
potential of an infinite straight wire along the z axis carrying current I. 

Solution 



We shall consider a long straight current along the z axis (—L < z < L ) and obtain the 
vector potential for cylindrical coordinate radii pCi, which is a good approximation for 
the case L oo. Then 


A(r) 


Mo 

47T 


Idr' 

~R 


/io/ f L dzz 

4?T yj £ + p 2 



ln(\/z 2 + p 2 + Z 4 


L 

-L 


(5.5) 

(5.6) 

(5.7) 
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£ ^/ ln /ygr t ± l \ 

4 tt \ vV /L 2 + p 2 -L/ 


(5.8) 


- IM)I , 
z —— In 
47T 


« z-In 


4?V 

ho-7 


Vi+p 2 /L 2 + r 
yi + P 2 /L 2 - ly 

1 + p 2 /2L 2 + 1 


= z : -In 

4tv \p 


1 + p 2 /2L 2 — 1 

(provided p <C L) 


(1st 2 terms of series expansion) 


2/2L 2 ; 


= z ^jr [ln(4L 2 ) - 2 ln(p)] 

A(r) = — z ln(p) + z C. (C is a constant.). 
27r 


(5.9) 


(5.10) 

(5.11) 

(5.12) 


Note that the value of C has no effect on the magnetic field, and so it is convenient to write 
A(r) = -z(/i 0 //27r)ln(p). 


5-3 We can add the gradient of a scalar field U (r) to A(r) without changing B(r). This is 
called a gauge transformation. The “gauge” of the vector potential is determined by the 
value of V • A. Show that in magnetostatics V • A = V 2 f/(r). 

Solution 


V • A(r) = V 

= ho /' 
47T J 




V • ( 1^2 ) dV + V 2 C/(r). 


(5.13) 

(5.14) 


We shall show that the integral is zero for a finite current distribution, i.e. for J(r —» oo) = 
0 . We can use the product rule V • (aV) = Va • V + aV • V to expand the integrand 


V 


J(v; 

R 


= V ^—) • J(r 7 ) + 0 (diff. is w.r.t. unprimed coords.) 
= -V' ( 1 ) • J(r') ( VR n = -V’R n = nR n ~ 1J R). 


(5.15) 

(5.16) 

(5.17) 
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Next use the same product rule but for primed coordinates V 7 • (aV) = V 7 a • V + aV 7 • V 
to get 


V 7 - 



= v '(s)' j(r ' ) + ^ v ' ,(r,) ' 

= V'(i).J(r' )+ 0, 


since charge conservation in magnetostatics requires V • J = 0. 
Then 


/ v ' 

Jail space 



d 3 r' = - [ V' • 

Jail space 



dV, 


s at 


J(Q 

R 


dS\ 


= 0 , 


since J(r 7 —>> oo) = 0 for a localised charge distribution. Hence, 
V • A(r) = V 2 U(r). 


(5.18) 

(5.19) 


(5.20) 

(5.21) 

(5.22) 


(5.23) 


In magnetostatics, it is convenient to choose U( r) such that V • A = 0 (Coulomb gauge). 


5-4 (a) Find the vector potential of the constant magnetic field Bq = (b) check 

that the vector potential you find does give the desired magnetic field, (c) find V • A and 
check it is what is expected in magnetostatics. [Hint: first find the vector potential A^(r) 
of the simpler constant field B® z by writing down the components of V x AW in Cartesian 
coordinates before appealing to the symmetry of the problem, and then integrating.] 

Solution 


(a) We first attempt to find A^(r) which must satisfy 

(d_A dAj z) \ = o 
dy ) 


(5.24) 
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A^)(r) must be symmetrical about the z axis which suggests that 

1 ft A^ 1 

dx 2 z ’ dy 2 

Integrating we get 

4 *’ = \ B ° 4 2> = -^»- 

Hence, 

AW = I(5°*y- J B«yx). 

Similarly, or by cyclic permutation of x, y and z above, 

A^ = ^(B° x yz-B° x zy), 

A (2/) = i (S° 2 ix- B°xz) , 


(5.25) 


(5.26) 


(5.27) 


(5.28) 

(5.29) 
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so that the vector potential of Bo = (B®, By,B®) is 

A = (aW + A(»4aW), (5.30) 

= \ [x ( B yZ ~ B ° z y ) + y ( B° z x - B° x z) + 2 (£°y - B° y x)] , (5.31) 

= ^B 0 x r = -ir x B 0 . (5.32) 

This A(r) is not unique - we could add the gradient of any scalar field to A(r) without 
changing B(r). 

(b) Taking curl of the vector potential found above, this time using index notation, 

V x A = -tv x (r x B°), (5.33) 

[V X A]j ^ £ ijk^ j £ klm'^lB m , (5.34) 

2 £ kij £ klmyjflB m , ( 5 . 35 ) 

= -^(SuSjm ~ S im Sji)VjnB^, (5.36) 

= -^uSjmVjriB^ - SimSjiVjnB^), (5.37) 

= - l -{V m riB° m -V l r l B*), (5.38) 

= -\{r^/ m B Q m + B° m V m ri - nVjfl? - 5?V,r,), (5.39) 

= -i(0 + £°-0-3£°), (5.40) 

= -^("2 Bf), (5.41) 

[V x A]* = Bl (5.42) 

.-. V x A = B° (as required). (5.43) 

We have used V • B° = 0, V m r^ = B° constant, and V • r = 3 used in Eq. 5.39 above. 
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(c) the divergence is 


-iv-(rxB°), 

(5.44) 


(5.45) 

+ 7j-Vi.Bg), 

(5.46) 

+ 0, 

(5.47) 

— 2 £ kij V iVj, 

(5.48) 

r] fc , 

(5.49) 

0 (because V x r = 0), 

(5.50) 


and the vector potential is seen to satisfy Coulomb gauge. 


5-5 A steady current / flows down a long cylindrical wire of radius b. Find the magnetic field 
both inside and outside the wire. 

Solution 



We define the z axis to correspond the axis of the wire, and point out of the screen/page. 
Then assuming the current density is constant inside the wire 



I /{rib 2 ) z 

0 


(0 < p < b) 
(p > b) 


(5.51) 
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From symmetry arguments B must be azimuthal, i.e. B(r) = B(p)<p. We apply Ampere’s 
law to loops ri and r 2 in the diagram, 


(J) B • dv — /^O^encl* 

(5.52) 

27rp 1 B(p 1 ) = (loop Ti), 

(5.53) 

2np2B(p2) = pol (loop r 2 ). 

(5.54) 

. B ( r \ _ [ {popl/^b 2 ) (j) (0 < p < b), 

1 (/ioT/27rp)0 (p>b). 

(5.55) 


5-6 A semi-infinite solenoid of radius a, has n turns per unit length, extends from z = —00 to 
z — 0 along the z axis and carries current I in the +</> direction. Magnetic flux is confined 
to the solenoid, but emerges isotropically from its end at the origin as shown below. 

(a) On the cone of half-angle 6 with apex at the origin there is a circular loop (as shown) 
with all points on the loop being at distance r a from the origin. Find the magnetic 
flux passing through this loop. 
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semi-infinite solenoid 



z 


(b) Find the magnetic vector potential at the point (r, 0, </>), and take it’s curl to find the 


magnetic field. [The expression for magnetic flux through a loop in terms of the vector 
potential may be useful here.] 

Solution 

a) The magnetic flux emerging from the end of the solenoid is the flux inside the solenoid, 
i.e. 


$0 = (tM) nl)(ira 2 ). 


(5.56) 


The flux emerges isotropically from the end of the solenoid at z = 0, so we shall need the 
solid angle subtended at the origin by the circular loop, which is 



‘9 


sin O'dO' = 2tt[— cos^Jq = cos#) — (—1)] 27r(l — cos #). 


(5.57) 


Then the magnetic flux through the circular loop is 



(5.58) 


(b) The magnetic flux through the loop is equal to the line-integral of the vector potential 
around the loop 



(5.59) 
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From symmetry arguments, since the current around the solenoid is only in the direction, 
the vector potential must also be in the <p direction, and since the loop has radius r sin # 


A^r , #, (j)){2nr sin #) 

■ A (r,e,(f>) 


= 2 ^ “ cos6 0 $ o, 

= 1 (1 -cosfl) 

2(2nrsm6) 0 

(1-COS0) 2 

= ~r, - (P 

(47r r sin 9) 


(5.60) 

(5.61) 

(5.62) 


Taking the curl, 


B = 


1 

r sin# 


' d_ 
89 


(sin 9 A'f) 


dA 0 1 _ 1 

“TT r + “ 
o<p J r 


1 dA r 
sin 9 d(j) 


d_ 

dr 


(r Ap) 


e 


+ 


i 

r 


d_ 

dr 


(r A e ) 



rsm9o9 r or 

1 d (1 — cos#) ^ 

=-^-4nr 4- 0 

r sin 6 dO (47rr) 

1 sin# T ^ 

= : 

r sin # (47rr) 

4>o ^ 

47rr 2 ’ 

fio(nlna 2 ) ^ 

4irr 2 


(5.63) 

(5.64) 

(5.65) 

(5.66) 

(5.67) 

(5.68) 


This field has similar form to the electric field of a point electric charge. Hence, the end 
of a semi-infinite solenoid appears as if it were a magnetic monopole of “magnetic charge” 
nlna 2 . 


5-7 By taking the curl of the vector potential for a magnetic dipole with moment m located 
at the origin, find it’s magnetic field using index notation. 

Solution 
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Using index notation, 


B(r) = gv X ( m x i) . 

p _ MO yy / ^kln^l r n\ 

Ul ~ UV £ ij k [ ^3 ) ’ 


47T 

Mo 


, ^kji^kln \ q 

47T J V H 


r " 

Vjmir n 


+ mir n Vj -3 ) , 


Mo 

^kji^kln 

Mo 


mi'Vjr, 


3' n 


r 3 +m,r n l-3^J -e,- 


mi 6 


— 47T ^ il ^ n ^ r 3 + m / r n( 3)rj^ 5 

Mo ^ mi5 nn TrijSji 3( TH/v n v rL TfijTiT j) 

47T 


Mo /3m* - m; 3[mjr 2 - (m • r) r*] 


47T \ H 
Mo /3(m • r)(r • e*) — m, 

-Di = T— 


B( r ) = Mo /^3(m • r)r — m 


(5.69) 

(5.70) 

(5.71) 

(5.72) 

(5.73) 

(5.74) 

(5.75) 

(5.76) 

(5.77) 
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5-8 A circular current loop in the xy plane has radius a and is centred on the origin. It carries 
current / in the (^-direction. There is a uniform magnetic field B(r) = £>o(cosdz + sindy) 
present. By integrating the torque dN = r' x dF on line element dr' of the current loop 
at r', find the torque on the entire current loop. Compare your result with the result you 
would get by first finding the current loop’s dipole moment, and then applying the formula 
for the torque on a magnetic dipole in a magnetic field. 

Solution 


Z A 


y 


A point on the loop and the corresponding line element on the loop are 



r' = a(cos 0 x + sin 0 y), 
dr' = ad(/)(— sin^x + cos^y). 

The force on that line element is 


dF = Idr' x B, 


= 1 


x y z 

—a d<j) sin 0 a d</> cos 0 0 


0 Bq sin 6 Bq cos 9 

= IaBo d0(cos0cosdx + sin^cosdy — sin (j) sin 9 z). 


(5.78) 

(5.79) 

(5.80) 

(5.81) 

(5.82) 
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The torque on the line element is 


x (dF), 



X 

y 

z 

acoscf) 

a sine/) 

0 

IaBo d<p cos (j) cos 6 

IaBo dc/) sin (j) cos 0 

—IaBo dc/) sin (j) sin 9 


= la 2 Bo d<j) [— sin 2 <fi sin 0x + cos (j) sin </> sin 6 y 

+ (cos (j) sin (j) cos 6 — sin </> cos (j) cos 9) z], 

= I a 2 Bo dc/) (— sin0 sin 2 </>x + cos (j) sin (j) sin 0 y). 


(5.83) 

(5.84) 


(5.85) 

(5.86) 


When integrating over <fi we note that cos (j) sin (f) is an odd function and its integral from 0 
to 27 t will be zero, and so 


N —la 2 Bo sin# x 



sin 2 cj) d(f) 


= —la 2 Bo sin# x 
= —l7ia 2 Bo sin# x. 



1 

- sm 
2 


cos 



Given that the dipole moment is m = In a 2 z, we expect a torque 


N = m x B, 

= —Ina 2 Bo sin# x 

as just found using the force law. 


(5.87) 

(5.88) 

(5.89) 


(5.90) 

(5.91) 
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6 Magnetism of materials 


6-1 Derive the following which are needed to obtain the magnetisation currents of a magnetised 
object: (a) Identity for V'i? -1 , 

V'RT 1 = +R~ 2 R. (6.1) 

(b) Product rule for V x (aF), 

V x (aF) = (Va) x F + aV x F. (6.2) 

(c) Corollary to Gauss’ Theorem, 

[ V x F d 3 r = - l F x dS. (6.3) 

Jv Js 
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Solution 


(a) First, for the identity for V'i? 1 we shall start by deriving the more general case, 

V'R n = [(x - x'f + {y- y') 2 + (z - z') 2 ] n/ 2 

+ y ^7 [(* - x ') 2 + ( y ~ y'f + ( z ~ z ') 2 ] n/2 

+ [(* - x'f + {y- y'f + (z - z ') 2 ] n/2 , (6.4) 

= x [(® - x'f + {y- y') 2 + (z- z'f ] (n/2_1) 2 (a; - x')(-l) 

+ y | [(* -z ') 2 + (y -y ') 2 + (*- X) 2 ] (n/2 - 1} 2( y - y '){-i) 

+ z | [(® - ®') 2 + (y - y ') 2 + (2 - ^') 2 ] (n/2_1) 2 (z - z')(-i) : (6.5) 

= ~n [x(x - x') + y(y - y') + z(z - z')] R n ~ 2 , (6.6) 

= -nKR n ~ 2 , (6.7) 

.-. V'R n = -nR n ~ x R. (6.8) 

Thus, for this exercise, 

V'R~ 1 =+R~ 2 R. (6.9) 

(b) The product rule for V x (aF) is derived as follows using index notation, 

[V x (aF)]* = E ijk S7j(aF k ), (6.10) 

= £ijk(F k \7ja + aVjFk ), ( 6 . 11 ) 

= &ijk(^F)jF k -\- a£ij k VjF k , (6.12) 

[V x (aF)]* = [(Va) x F]* + [aV x F]*. (6.13) 

.-. V x (aF) = (Va) x F + aV x F. (6.14) 
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Finally, re-arrange the product rule: 


— (Vo) x F = aV x F 

- V x (oF), 

(6.15) 

F x (Va) = aV x F 

- V x (oF). 

(6.16) 


(c) We shall prove the corollary to Gauss’ Theorem by applying Gauss’ theorem for the 
vector field c x F(r) where c is a constant vector is 



(c x F) • dS. 


(6.17) 


Then we rearrange the right-hand side using the scalar triple product rule (c x F) • dS = 
(F x dS ) • c, and on the left-hand side we can integrate by parts using the product rule 
V • (c x F) = (V x c) • F — (V x F) • c, giving 


[ [(V x c) - F- (V x F) -c] d 3 r = l(Fx dS) • c, 

Jv Js 


[ V x F d 3 r 

• c = 

® F x dS 

Jv 


Js 


C, 


I V x F dr = -| Fx dS, 


as V x c = 0 for c being a constant vector. 


(6.18) 

(6.19) 

( 6 . 20 ) 


6-2 A thin disc of magnetised material is coincident with the xy plane. It is of thickness s and 
radius a and has magnetisation M = Mqz. Find the magnetisation current, and from this 
find the magnetic dipole moment of the disc. Compare this with what you would get by 
multiplying the disc’s volume by M. 

Solution 
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The surface magnetisation current is 

K mag = Mxn = M 0 zxp = M o 0. (6.21) 

The net magnetisation current around the disc’s circumference is / mag = sif mag = sMq, 
and so the dipole moment is 

m (7ra 2 )/ mag z 7ra 2 sM. (6.22) 

This is just the volume multiplied by the magnetisation field. 

6-3 Consider a permanent magnet in the form of a short cylinder of radius a extending along 
the z axis from z — —L to z — -\-L and having uniform magnetisation M = Mq z. (a) Find 
B and H at all points (0,0, z) on the cylinder’s axis, and plot B(0, 0, z) and H(0,0,z) vs. 
z. (b) Discuss whether the result obtained in part (a) obeys Ampere’s law for H in integral 
form. 
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Solution 



(a) Because the magnetisation field is uniform, there is no magnetisation volume current. 
The magnetisation surface current K mag (r) = M(r) x n is zero on the two ends, and 
K mag (r) = Mq<P on the cylindrical surface. 

The Biot-Savart law gives the magnetic field due to an arbitrary surface current distribu¬ 
tion: 


B(r) = [ [Kmae(r,) f] X ^ (6.23) 

47 T J R 2 

The diagram shows the contribution AB to the magnetic field at (0,0,2:) due to the surface 
magnetisation current in a small patch of the surface making up part of the strip of thickness 
dz' at z'. The components of AB due to the surface currents in different patches around 
the strip which are not in the 2 : direction will cancel each other out. This leaves only a 
2 :-component for the contribution d'. B, due to the entire strip, so that 


g?B(0, 0, z) 


di. B(0,0, z) 


Mo (Kmdigdz )(2na) ^ 

— COS u Z 

47t R 2 

Mo ( Mdz r ){2na ) a ^ 

4 ^ E? R Z ’ 


fi^McPdz' _ 
2[0-z') 2 + a 2 ] 3 /2 Z ' 


(6.24) 

(6.25) 

(6.26) 


71 


Download free eBooks at bookboon.com 



















Essential Electromagnetism: Solutions 


Magnetism of materials 


Hence, integrating over the entire cylindrical surface, 

fi^Ma^dz 1 


B(0,0, z) 


-L 


_ L 2[(z — z') 2 + a 2 ] 3 / 2 


z, 


B(0,0, z) = 


/ioM 

2 

Mo M 


z + z' 


n L 


y/(z + z') 2 + a 2 


J —L 


z -\- L 


z — L 


\/(z + L) 2 + a 2 \J{z — L) 2 + a 2 


z. 


Now, H = B//i() — M, and so 


H(0,0,z) = < 


M 0 

z~\~L 

z—L 

2 

< 

_\J (z+L) 2 +a 2 

y/(z-L) 2 +a 2 _ 

M 0 

z-\-L 

z—L 

2 

s. 

yO+L) 2 +a 2 

yj(z—L) 2 -\-a 2 _ 


z (\z\ > L ) 

z — Mqz (|z| < L) 


B and H are plotted below. 



z/L 


(b) Ampere’s law in integral form is 

H • dx If end. 


(6.27) 

(6.28) 

(6.29) 


(6.30) 


(6.31) 
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where If, e nc i. is the net free current through loop I\ Since there are no free currents the 
integral must be zero for any closed path. We only know H on the z axis, and at |r| = oo 
where H must be zero because the source of magnetic field, in this case the magnet, is 
localised near the origin. But we can construct a closed loop which has as part of it the 
entire z-axis as follows: (0, 0, — oo) to (0, 0, +oo) to (0, +oo, +oo) to (0, +oo, —oo) and back 
to (0,0, — oo). The integrand is zero except along the z-axis, so in this case Ampere’s law 
in integral form is satisfied provided 

/ oo 

H z {0,0,z)dz = 0. (6.32) 

-OO 

Examining the plot of H z vs. z above, it appears that the integral is indeed zero, as could 
easily be checked by numerical integration. 
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6-4 A cylindrical rod of radius a and length h a is permanently magnetised along its length 
which coincides with the z direction, i.e. M = Mqz. (a) Find the surface magnetisation 
current K mag (r), and use it together with Ampere’s law to find B and H inside and outside 
the rod (assume h —oo). (b) The rod is now bent into a circle of circumference (h + 2 L) 
such that there is an air gap of width 2 L < a. Plot B and H along the axis of the magnet 
in the vicinity of the air gap for the case of a — 0.5 and L = 0.2. 

Solution 



(a) The surface magnetisation current density is K mag (r) = M(r) x n, hence 

K mag (a, 0, z) = M 0 zxp = M o 0. (6.33) 

This surface magnetisation current is similar to the current in a tightly-wound solenoid, 
and the magnetic field inside the rod can be calculated in the same way using Ampere’s 
law j> B • dr = /ioAmci- From the symmetry of the problem B inside the rod can only be in 
the z direction. 

For Amperian rectangular loop T (see diagram) with one side of length Sz inside the rod 
at cylindrical radius pi and one outside at cylindrical radius p 2 


[B z (pi,(/),z) - B z (p 2 ,(/),z)\ Sz = pqMqSz. (6.34) 

That this is independent of p 2 and applies equally to p 2 —> oo (where B = 0) tells us that 
B = 0 outside the rod. Again, since the integral is independent of p\ the magnetic field 
inside the rod is constant, 


B(p<a,0,z) = /hqMqz = p 0 M. 


(6.35) 


74 


Download free eBooks at bookboon.com 









Essential Electromagnetism: Solutions 


Magnetism of materials 


Finally, 

H = ( — j - M = 0 (6.36) 

V/W 

everywhere. Note that this result is for an (unrealistic) infinite magnetised rod, and that 
near the two ends of the rod B would be different, and H would be non-zero. 

(b) To find B and H on the axis of the magnet near the air gap we can use the information 
that h a and assume the magnet in this region is approximately straight, with its axis 
being along the z- axis, and with the air gap extending from z = — L to z = +L. In that 
case we can imagine that a short magnet of length 2 L has been removed from an infinite 
magnetised rod. 

Using the principle of superposition, we can get B in the vicinity of the air gap by subtract¬ 
ing the field of the short magnet of length 2 L from the field of an infinite straight magnetised 
rod with no air gap. For the short magnet we can use the results from Exercise 6-3 for the 
magnetic field of the short magnet. Thus, 


B(0,0,z) = [jqMqz - 


MoM 0 


z L 


z-L 


y/(z + L) 2 + a 2 y/(z - L) 2 + a 2 


z. 


(6.37) 


Now, H = B/^o — M, and so 


H(0, 0, z) = < 


Mo 


z-\-L 


z—L 


2 yj (z+L) 2 +cl 2 yj ( z-L) 2 +cl 2 

z-\-L z—L 


M 0 z - ^ 


V (*+£) 2 +« 2 


(\z\>L) 
z (\z\ < L ) 


(6.38) 


B and H are plotted below. Note that the area under the plot of H z ( 0, 0, z) vs. z appears 
to be zero in agreement with Ampere’s law in integral form for the case of no free currents. 
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6-5 Consider a permanently magnetised sphere of radius a with uniform magnetisation M(r) = 
Mqz. (a) Find the surface magnetisation current density, and use this to find the magnetic 
dipole moment of the sphere. Compare this with what you expect given the volume of the 
sphere and the magnetisation field, (b) Find B and H at the centre of the sphere. 

Solution 



(a) The surface magnetisation current density is 


K mag (a, 0, z) — Mxn = Mqz xr = Mq sinO (j). 


(6.39) 


The magnetic dipole moment of the surface magnetisation current distribution is 

m = i f r x K mag (r )dS. (6.40) 

From symmetry arguments, the dipole moment must be m = m z z where 


If 1 — 

z - / ar x (MosinO <p)27ra 2 d(cos6) 

2 J- 1 

(6.41) 

2tvo^Mq f 1 ^ 

--- / sin^z • (— 9) d(cos6) 

2 J- 1 
f 1 

(6.42) 

7ra 3 Mo J sin 2 6 d(cos 6) 

(6.43) 

7ra 3 Mo J (1 — cos 2 6) d( cos 6) 

(6.44) 

4 Q 
-7ra 3 M 0 

(6.45) 


as expected. 
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(b) We can use the Biot-Savart law to obtain the magnetic field at the centre of the sphere 

B(0,0,0) = ^ <f Kmas( p r ? X ~ dS'. (6.46) 

4jV J R z 

From symmetry, this must be B(0, 0,0) = B( 0, 0, 0) z where 


5(0,0,0) = 

- Mo f 1 (M 0 sin00)x(-r) 2 

z-— / - tz - 27rad(cos(J), 

47r J_ 1 a z 

(6.47) 

= 

/ sin#z • (— 0) d(cos6), 

2 J- 1 

(6.48) 

= 

MoM 0 f 1 . 2 Q 1, m 
—-— / sin ^(cos#), 

2 J- 1 

(6.49) 

= 

° f (1 — cos 2 6) d(cos$), 

2 J- 1 

(6.50) 

5(0,0,0) = 

2 

-/IqMq. 

(6.51) 

B(0,0,0) = 

2 

(6.52) 


and is in the same direction as M. 

Now H = B//io — M, so that H(0, 0, 0) = — and it is in the opposite direction to M. 


6-6 Consider the hysteresis loops of the magnetically-soft iron-based amorphous alloy and the 
magnetically-hard alloy of iron, aluminium, nickel and cobalt shown in Chapter 6. 

(a) Estimate the work done to bring 1 cm 3 of each material through one cycle of the 
hysteresis loop, (b) Two transformers operating at 50 Hz have magnetic cores of volume 
100 cm 3 (one of each type of material) and are (unwisely) operated at a current at which 
saturation occurs. How much power is lost as heat in each case? [You could print the 
hysteresis plots and estimate the area by drawing over it a grid and measuring by hand 
sufficient points on the graph to get within say 20% accuracy for the area.] 

Solution 
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-400 -200 0 200 400 

H (kA m' 1 ) 



-400 -200 0 200 400 

H (kA m' 1 ) 


Grid lines have been drawn over the hysteresis plots above. The area of the upper half of 
the hysteresis loop (for positive %) is identical to that of the lower half. Read off the 
(horizontal) “widths” in H at “heights” Bm = 0, 0.1, 0.2, ... (T). The approximate values 
are tabulated below. 



American online 
LIGS University 

is currently enrolling in the 
Interactive Online BBA, MBA, MSc, 
DBA and PhD programs: 


► enroll by September 30th, 2014 and 

► save up to 16% on the tuition! 

► pay in 10 installments/2 years 

► Interactive Online education 

► visit www.liqsuniversitv.com to 
find out more! 


Note: LIGS University is not accredited by 
nationally recognized accrediting agency 
by the Ub Secretary of Education. 

More info here. 


79 


^3 


Click on the ad to read more 


Download free eBooks at bookboon.com 






























































































































































































Essential Electromagnetism: Solutions 


Magnetism of materials 


Bm (T) 

iron-based alloy 
AH (kA m -1 ) 

Alnico 
AH (kA m -1 ) 

0.0 

30 

120 

0.1 

25 

120 

0.2 

25 

120 

0.3 

20 

120 

0.4 

20 

120 

0.5 

15 

120 

0.6 

15 

120 

0.7 

15 

120 

0.8 

10 

120 

0.9 

10 

120 

1.0 

0 

120 

1.1 


115 

1.2 


100 

1.3 


90 

1.4 


0 

sum: 

185 

1625 


The areas of the two loops are approximately: 


iron-based alloy: j) Bm^H = 2x (185 kA m 1 ) x (0.1 T) = 3.7 x 10 4 J, (6.53) 
Alnico: £ B M dH = 2 x (1625 kA m" 1 ) x (0.1 T) = 3.2 x 10 5 J. (6.54) 

Now, this is for a sample volume of 1 m 3 . For a magnetic core of volume 100 cm 3 =10 -4 m 3 , 
the energy to take the sample around one cycle is then 3.7 J (iron-based alloy) or 32 J 
(Alnico). 

If the sample is used in a transformer operating at 50 Hz with the magnetic field saturating, 
the energy lost to heat in one second is 50 times the energy for one cycle, i.e. 185 W (iron- 
based alloy) or 1.6 kW (Alnico). If the core were made of Alnico, the rate of heating would 
be similar to that of a domestic electric room heater! 
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